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ABSTRACT 


This  note  considers  the  problem  of  detecting  the  presence  of  pseudonoise 
waveforms.  The  detector  is  assumed  to  have  no  knowledge  of  the  sequence  pro¬ 
viding  the  pseudorandom  phase  shift  keying.  It  is  found  that  the  structure 
and  performance  of  the  optimum  detector  is  nearly  the  same  as  a  radiometer 
detector  with  proper  prefiltering.  In  the  special,  but  important,  case  when 
the  received  phase  of  a  binary  pseudonoise  signal  is  constant  over  the  de¬ 
tection  interval  it  is  found  that  the  optimum  detector  performance  is  3  dB 
better  than  the  radiometric  detector. 
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I.  INTRODUCTION  AND  CONCLUSIONS 


In  this  note  the  structure  of  the  optimum  detector  of  pseudonoise  wave¬ 
forms  is  derived  for  three  channel  models.  The  three  models  are: 

A.  Phase  known  to  the  detector, 

B.  Phase  unknown  to  the  detector  but  constant  over  the  detection 
interval , 

and  C.  Phase  random  from  chip  to  chip. 

In  each  case  the  performance  is  analyzed  fcr  the  practical  situation  of  a  low 
SNR  per  chip. 

The  received  waveform,  r(t),  is  observed  for  N  chips,  each  of  duration 

Tc.  Each  chip  waveform  is  of  the  form  ±  /  cos  (wt  +  0^  where  0,  is  the 

c 

phase  on  the  i-th  chip.  The  sign  is  unknown  to  the  detector.  In  addition  it 

is  assumed  that  additive  white  Gaussian  noise  of  single-sided  density  Nq  is 

present.  Statistics  sufficient  for  the  detector  to  make  a  decision  consist 

of  the  resultsjjf  correlating  each  chip  against  in-phase  and  quadrature 

references,  J ^  ^  eos(u)t)  and  J -  —  sin(wt).  These  will  result  in  two 
o  c  o  c 

statistics  for  each  chip,  r^  and  r^  normalized  to  unity  variance.  Note  that 

this  correlation  is  equivalent  to  passing  the  received  waveform  through  filters 
sin  (!*f/T  ) 

with  .  j j ^  amplitude  characteristics  but  90°  offset  phase  characteristics 

c 

and  sampling  the  outputs  every  Tc  seconds.  While  synchronization  is  assumed 
in  each  case,  this  is  not  a  severe  restriction.  Clearly  several  time  origins 
may  be  examined,  one  of  which  will  be  close  to  the  results  presented  here. 


The  false  alarm  rate  would  consequently  rise  -  but  an  insignificant  gain  in 
required  detector  SNR  will  compensate  for  this. 


1 


The  behavior  of  the  optimum  detector  for  binar,  pseudonuiy,-  waveforms 
with  small  Ec/NQ  is  found  to  be  similar  to  a  radiometric  detector  preceded  by 
a  signal  spectrum  matching  filter.  However,  in  cases  of  practical  interest 
there  is  a  3-dB  improvement  in  SNR  compared  to  the  radiometer.  If  the  phase 
of  the  received  signal  is  known,  the  detector  SNR  will  be: 


2T 


P  2 
<ir> 


(1/T  )  VN 
c  o 

where  the  total  observation  time  is  T=NT  ,  and  the  received  power  is  P  *E  /T  . 

c  r  c  c 

The  false  alarm  probability  and  detection  probabilities  being  given  by  the 
following  trade-off  equations: 

Pf  -  Q(Y) 

P,  ”  Q(Y-d) 


*  2 
where  y  is  a  threshold.  If  the  phase  is  random  from  chip-to-chip  d  is 

reduced  by  3  dB.  (A  change  in  required  (Pr /Nq)  of  only  1.5  dB.)  The  case  of 

unknown  but  constant  phase  (which  clearly  must  be  between  the  two  previous 
2 

cases)  has  a  d  essentially  equal  to  the  known-phase  case. 

It  should  be  noted  that  higher  order  phase  modulation  (e.g.,  quadraphase) 
would  be  bounced  in  behavior  by  the  known  and  random  from  chip-to-chip  cases 


above . 


II.  KNOWN-PHASE  MODEL 

If  0^  is  known,  the  r^  terms  can  be  assumed  to  be  in-phace  with  the 
received  waveform,  r^  can  be  ignored.  The  likelihood  ratio,  A(r),  where  r 

*  2TT 

A  .  e-x  '2 

q(x>  =  /  — —  dx 

x  /2tt 


2 


refers  to  the  entire  set  of  r^  and  r^  observed  can  be  written  as: 


Mr)  *“JT 

N 


1  i  r-<ru->>/ 

T  -  e 

2  ffi  L 


-m)  /2  -(rj^+m)' 


-r2  /2 
1  rIi/z 
—  e 


where  m  *  / .  The  product  is  taken  over  the  N  chips  under  observation, 
o 

Each  numerator  term  is  just  the  probability  density  function  (pdf)  of  r^  which 
is  Gaussian  with  unit  variance  and  mean  ±m  with  equal  probability.  Canceling 
common  terms  and  factoring  terms  out  of  the  product  which  are  independent  of  r 
gives  the  simp.1  form: 


const 


tt[, 

N  L 


'  2rn"  ^  -2rii” 
e  +e 


At  this  point,  A(r)  can  be  manipulated  in  two  ways.  If  the  product  were 

„  2m£±rT. 


expanded  into  a  sum,  there  would  be  2  terms  each  of  the  form  e 


where 


the  sum  in  the  exponent  would  vary  over  all  the  2n  ways  that  signs  can  be 

chosen  to  combine  the  r^  terms.  This  shows  that  the  likelihood  ratio  for  this 

composite  hypothesis  testing  problem  is  the  sum  of  tb'’  likelihood  ratios  for 

N 

each  of  the  possible  2  received  waveforms,  a  well-established  fact  but  not  a 
practical  way  ot  designing  a  detector. 

A  more  interesting  manipulation  is  the  following: 

A(r)  =  const  |]  cosh  (2rI.m) 

N 

after  taking  the  logarithm  (neglecting  constants): 


L(r)  =  Y.  £n  cosh  (2r  .m) 
N  1 


3 


v-!  r»  j  ly  ^ -.t 


which  shows  that  2  filters  need  not  be  built.  Indeed*  L(r)  consists  of  a  non¬ 
linear  weighting  of  the  matched  filter  outputs  which  are  then  summed.  The  sum, 
L(r),  would  then  of  course  be  compared  to  a  threshold. 

K 

»*  / » 


The  most  interesting  practical  case  is  when  m,  /  zr^-  ,  is  small.  In 

x2  ° 

this  case,  the  approximation  An  cosh  (x)  «  —  can  be  used  to  show  that  (except 
for  constants): 

L(r)  a  T  r2  . 

N 

Hence  the  optimum  detector  is  a  square-law  or  energy  detector.  However  'half' 
of  the  received  energy  is  not  included  since  it  is  known  to  be  noise.  A 
structure  for  the  optimum  receiver  is  shown  in  Fig.  1. 


The  impulse  response  of  the  filter  shown  is: 

4 


h  (t) 

c 


N  T 
o  c 


COSuJt 


0  <  t  <  T 


0  elsewhere. 

The  performance  of  this  detector  follows  readily  from  noting  that  L(r)  is 
the  sum  of  many  independent  variables  and  hence  can  be  considered  Gaussian. 
Thus  all  that  is  required  for  performance  evaluation  are  the  means  and 
variances  under  noise-alone  (hypothesis  Hq)  and  signal-plus-noise  (hypothesis 
Hj)  conditions.  For  noise  alone: 

E  [L(r>]  =  N  E  [r2J 

=  N 


and 


4 


Fig.  1.  Optimum  detector-known  phase. 


*  THRESHOLD 
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Var  [L(r)] 


N  V.r  [r^] 


-E2[r2]) 


N  (3-1) 


which  all  follows  from  the  fact  that  r^  is  Gaussian,  zero  mean  and  normalized 
to  unit  variance.  Assuming  signal-plus  noise: 

E  [L(?)]  »  NE  [r2  3 


E  [L(?)]  »  NE  [r2J 

■  N  (1  +  m2) 
Var  L(r)  *  N  Var  [r2  ] 


-  N  (E  [rJJ  -  E2  [r2^  ) 

-  N  (m4  +  6m2  +  3  -  (1  +  m2)2) 

=•  N  (4m2  +  2) 

which  follow  from  the  fact  that  r^  is  a  zero  mean,  Gaussian  with  unit 

2 

variance  added  to  ±m.  Since  moments  of  r^  are  being  computed  the  sign  of  m 
is  irrelevant . 

Hence  when  signal  is  present  the  mean  of  L(r)  increases  by  Nm2,  the 

2 

variance  remains  2N  when  terms  of  order  m  are  neglected.  The  "signal-to- 

o 

noise  ratio"  of  this  detector,  d  ,  is  then: 

2 

cj 2  =  (A  mean) 
variance 


F  2 

-  2"  fch  ■ 
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Relating  this  to  the  received  power,  P  -  E  /T  and  the  total  observation  tine, 

r  c  c 

T  *  NTc  gives: 


2T 


P  2 

<5T> 


(1/Tc)  'K0 

which  is  3  dB  better  than  a  radiometer  of  integration  time  T  trying  to  detect 
a  signal  in  bandwidth  W  -  1/Tc .  The  3-dB  gain  comes  about  because  of  the  known 
phase  (as  will  be  shown  below).  It  shculd  be  noted  that  this  detector  actually 
uses  the  full  signal  energy  by  correlating  with  chip  reference  waveforms.  A 
detector  which  consisted  merely  of  a  rectangular  filter  of  bandwidth  1/T  H? 

followed  by  a  sampler,  squarer  and  adder  would  have  about  2  dB  poorer  performance 

2 

(i.e.,  d  reduced  by  2  dB)  since  1  dB  of  the  signal  power  would  be  filtered  out. 


III.  UNKNOWN,  CONSTANT  PHASE 


The  case  when  the  0^  are  all  equal  (to  0),  but  unknown  proceeds  in  a 
similar  way.  The  likelihood  ratio  can  be  written  as: 


.  -(r,.-m  cos  0)2/2-(r  .-m  sin  0)2/2 

—  re  11 

4tt  le 


*<?>  -  TT 

N 


+  C0S  ®)2/2“(rQi‘Hn  9in  0)2/2n 

+e  ^  __ _ _ _ J_ 


1  ”(rIi  +  rQi  )/2 

•  ■**  p 

2  it  e 


const 


m(rn  cos  0+r  sin  0)  -nKr^  cos  e+rQi  sin  ^ 
e  +e 


=  const  |P  cosh  m(rj^  cos  9  +  rn-t  s*n 


Qi 


where  the  whole  expression  must  be  averaged  over  P  (which  is  assumed  to  have  a 
uniform  distribution). 
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If  the  product  is  expanded  into  a  sum  before  averaging,  2N  terms  will 
appear  of  the  form: 

m  [(E  ±  r^)  cos  0  +  (E  ±  r^)  sin  0]° 

where  again  the  E+  notation  refers  to  all  combinations  of  sign  used  in  the 
audition.  Then,  using  the  Bessel  relationship, 

1  /27r  a  cos  0  +  b  sin  0  I  a2+b^  ) 

-rr  /  e  d0  o 

2TT  *Q 

performing  the  average  over  0  will  lead  to  a  sum  of  2N  terms  of  the  form 
I^mX^fj^)2  +  (E+r^)2)  .  This  again  is  the  sun  of  the  likelihood  ratios 

for  each  possible  2N  waveforms  after  envelope  detection.  (The  Bessel  function 
argument  is  the  envelope.) 

The  more  interesting  approach  for  small  m  is  again  to  use  the  expansion 
2 

for  in  cosh  x  w  y-  to  show  that 

Xt,r)  *  exp  [I  in  (cosh  (mCr^  cos  0  +  r^  sin  0)))JD 

- 2 - 6 

«  exp  y-  [Efr^  cos  0  +  r^  sin  0)2] 

2 

*  exp  y-  [(Ir^Jcos2  0  +  (Er^Jsin2  0  +  2(ErIirpi)cos  0  sin  0] 

- - - e 

*  exp  y-  (A  cos2  0  +  B  sin2  0  +  2  C  cos  0  sin  0) 

where  A,  B,  C  represent  the  three  sums.  The  average  over  0  can  ue  performed 
by  integrating  over  (0,  2n)  which,  after  some  trigonometry  yields: 


8 


where 


A(?) 


m2  (A+B)  / 4  /-2lT  m2/2  [(A-B)/2  cos  28  +  C  sin  29]  JQ 

2tijq  d0 

em2(A+B)/4  /SZ77) 


A 

B 


£  r 
£  r 


2 

Ii 

2 

Qi 


C  *  £  r 


IirQi* 


Note  that  the  only  approximation  used  to  derive  the  above  is  chat  of  small  m. 

If  it  is  now  assumed  that  the  argument  of  the  Bessel  function  is  large  (greater 

than  «3)  its  asymptotic  expansion,  I  (x) - —  ex,  can  be  used.  Ignoring  the 

A  term  and  other  constants,  then  taking  the  logarithm  of  X(r)  yields 

i  • 

uh  «  -aai  +  c2  . 

.2..  .4  , 

It  will  be  presently  shown  that  the  approximation  of  the  large  argument  is  a 
good  one.  In  any  case,  the  above  expression  defines  a  good  receiver  structure 
even  if  not  the  optimum  one.  After  some  manipulation  the  receiver  can  be 
shown  to  be  of  the  form  shown  in  Fig.  2.  The  detector  has  an  interesting 
interpretation.  The  top  branch  is  essentially  a  radiometer  preceded  by  a 
filter  matched  to  the  signal  spectrum.  The  output  of  the  squarer  also  contains 
double -frequency  energy.  The  bottom,  double  branch,  is  used  to  detect  double - 
frequency  energy  coherent  over  the  measuring  interval.  It  is  similar  to  the 
structure  of  a  squaring  loop  which  utilizes  the  fact  that  squaring  a  biphase 
signal  plus  noise  strips  off  the  biphase  modulation.  As  will  be  shown  below, 


9 


the  low-pass  and  double-frequency  noise  terms  are  independent,  thus  it  is  not 
surprising  that  both  bands  should  be  used  for  detection. 

The  analysis  of  this  detector  is  unfortunately  rather  tedious.  Hence 
only  a  sketch  of  the  derivation  will  be  given.  The  assumption  that  will  be 
made  is  that  A,  8  and  C  are  jointly  Gaussian  since  each  is  the  sum  of  many 
independent  terms.  The  following  second-order  statistics  can  be  derived  in  a 
straight-forward  manner  (assuming  hypothesis  Hj  and  mx  =  m  cos  0,  mg  =  m  sin  ti) 


E[A] 

-  N(l4®j) 

Var  [a] 

-  N(4mj  +  2) 

E  [B] 

-  N(l-hn^) 

Var  [b] 

-  N(4mp  +  2) 

E  [C] 

■  Nm^m^ 

Var  [C] 

■  N(l+mj  +  m^) 

The  terms  appearing  in  L(r)  are  A+B,  A-B  and  C  which  will  also  be  Gaussian. 

It  can  be  shown  that  these  three  variables  are  uncorrelated  (hence  independent  1 
and  that  Var  [a+b3  ■  Var  [A-b]  *  N(4+4m^).  The  mutual  lack  of  correlation 

shows  that  the  double-frequency  noise  terms  are  independent  of  the  low- 
frequency  terms  in  Fig.  2. 

L(r)  is  now  seen  to  be  the  sum  of  two  independent  terms,  which  is 


Gaussian,  with  mean  N  and  variance  N,  and 


/U 


+  ^,2  has  a  Raleigli 

4 


distribution  with  mean  /  N  under  H  .  Under  H.  the  second  term  has  a 

2  o  l 


Rician  distribution.  By  a  direct  calculation  it  can  be  shown  that  the  p.d.f. 


for  L(r)  assuming  Hq  is: 


Note  that  the  Bessel  function  argument  has  mean 
enough  to  justify  the  large-argument  approximation. 


J TTNn/ 


which  is  large 


11 


P  (£)  ,  JL  [~ (&-N)Q (  (N-A) /  Jj)  e-a-N)2/4N+  _e_ 

L  *je  L  7zz  rr. 


-(A-n)72n 


]• 


L  2vfi  L  /2N  /2ttN 

Since  one  i3  most  always  interested  in  low  false  alarm  probabilities,  it  is  the 

tail  of  this  distribution  which  is  of  most  importance.  This  is  the  region 

1  -x^/2 

where  H  »  N  and  the  approximation  Q(x)  ■  1  -  ^  e  for  x  «  0  can  be  used. 
Keeping  only  the  largest  term  gives: 

pm  «  ShSX  e-(-»)2/ 4N 

L  2/2N 

trotn  which  it  follows  that  the  false  alarm  probability  is: 


Pr(L>i.) 


1  -(&-N)^/4N 

- e 

»^N 


Thus  the  false  alarm  probability  behaves  as  if  the  test  statistic,  L,  was 
Gaussian  with  mean  N  and  variance  2N. 

^2 

Under  Hj,  the  mean  of  — ~  increases  by  -j-  while  the  variance  stays 

essentially  unchanged  at  N.  The  square -root  term  will  be  essentially  Gaussian 
2  ^ 

with  mean  — !•—  and  variance  N  as  long  as  »  1  which,  it  will  be  shown,  is 

the  condition  for  good  signal  detectability.  Thus  under  Hj,  L  is  again 

2  2 

Gaussian  with  variance  2N  with  the  new  mean  N  +  Nm  .  It  follows  that  d  can 

be  defined  in  the  same  way  as  in  the  known-phase  case  and  will  be: 

.2  Nm4 
d  -  — 


or 


2T 


Pr  2  2 

(t r)  d*  »  l. 


(1/T  )  VN 
c  o 

This  result;  i3  identical  to  that  for  known  9  which  is  reasonable  in  retrospect 
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v  JT^-TSW-* 


n^Trr^?\*^Trflf*7*If*i  ?-V’ 


since  it  can  be  shown  that  the  received  phase  can  be  estimated  with  standard 


deviation 


«  1 


radians. 


Hence  the  phase  estimate  will  be  quite  accurate  for  large  d  making  this  case 
equivalent  to  the  case  of  known  0.  (The  phase  estimate  with  this  accuracy 


is  given  by: 


I  2C 

*2  arctan  ^  .) 


IV.  PHASE  RANDOM  FROM  chip-to-chip 

In  this  case  the  0.^  are  assumed  to  be  uniformly  distributed  from  0  to  2tt 
and  independent  from  chip-to-chip.  The  likelihood  ratio  is  similar  to  that  of 
the  previous  case,  except  here  the  averaging  over  0^  is  done  for  each  term  in 
the  product  before  multiplying: 

A(r)  *  const  TT  cosh  inM..  cos  0 .  +  r  .  sin  ©i)  6i 
N  V 


*  const  7[  ^(m/rji  +  r^)  . 

N 

In  this  case  it  is  seen  that  there  is  no  structure  which  can  take  advantage  of 
the  2^  possible  transmitted  waveforms  since  they  are  completely  scrambled  by 
the  random  phase  shifts.  Taking  logarithms  and  neglecting  constants: 

L(r)  =  l  £n  +  rq7  • 

N 
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For  small  m  the  sufficient  statistic  can  be  simply: 


Hq,  noise  alone:  E  [L(r)]  -  2N 

Var  [L(r)]  -  4N 

Hj,  signal  plus  noise:  E  [L(r)]  *  N(2+or) 

Var  [L(r)]  «  N(4  +  4m2). 

Hence  for  this  case: 

2  Nm4 

d  "  T 

_  P  2 

-  -JL  (-£) 

l/T  'N  ' 
c  o 

which  is  essentially  the  same  ae  a  radiometer  detector  with  proper  pre-filtering 

2 

and  only  3  dB  worse  in  terms  of  d  than  the  known-phase  case.  (The  difference 

in  terms  of  P  /N  will  be  only  1.5  dB.) 
r  o 
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